Microscopic origin of diagonal stripe phases in doped nickelates 

Marcin Raczkowski 

Marian Smoluchowski Institute of Physics, Jagellonian University, Reymonta 4, PL-30059 Krakow, Poland ana 

Lahoratoire CRISMAT, UMR 6508 CNRS-ENSICAEN, 6, 

Boulevard du Marechal Juin, 14050 CAEN Cedex, France 



o 
o 

(N 

u 

< 



u 



o 



> 

oo 
en 

O 

o 



o 



X 



Raymond Fresard 
Lahoratoire CRISMAT, UMR 6508 CNRS-ENSICAEN, 6, 
Boulevard du Marechal Juin, 14050 CAEN Cedex, France 

Andrzej M. Oles 
Marian Smoluchowski Institute of Physics, Jagellonian University, Reymonta 4, PL-30059 Krakow, Poland 

(Dated: February 6, 2008) 

We investigate the electron density distribution and tlie stability of stripe phases in the realistic 
two-band model with hopping elements between Cg orbitals at Ni sites on the square lattice, and 
compare these results with those obtained for the doubly degenerate Hubbard model with two 
equivalent orbitals and diagonal hopping. For both models we determine the stability regions of 
filled and half-filled stripe phases for increasing hole doping x = 2 — n in the range of a:: < 0.4, using 
Hartree-Fock approximation for large clusters. In the parameter range relevant to the nickelates, 
we obtain the most stable diagonal stripe structures with filling of nearly one hole per atom, as 
observed experimentally. In contrast, for the doubly degenerate Hubbard model the most stable 
stripes are somewhat reminiscent of the cuprates, with half-filled atoms at the domain wall sites. 
This difi'erence elucidates the crucial role of the off-diagonal Cg hopping terms for the stripe formation 
in La2-iSra;Ni04. The influence of crystal field is discussed as well. 

PACS numbers: 75.30.Kz, 71.10.Fd, 75.10.Lp, 75.50.Ee 



I. STRIPE PHASES IN NICKELATES 



Stripe phases are one of the most exciting phe- 
nomena of modern condensed matter physics. They 
have been observed in a variety of systems, includ- 
ing nickelatesj^i^i'^i'^i^i^/i^ cuprates;Qd°i"d^i^'^d'^i^^ 
and manganites^SiiLiS Among them, layered 
La2-:,Sr:,Cu04 (LSCO), La2-^-yNdySra;Cu04 (Nd- 
LSCO), La2-a;Sr^Ni04 (LSNO), and La2Ni04+5 (LNO) 
compounds play plausibly the most prominent role. 
However the similarity between them is superficial only, 
and the stripes in the cuprates differ from the stripes 
in the nickelates in many respects. For instance, they 
are dynamical in the former, and static in the latter. 
In addition, in Nd-LSCO^'^Siii-iSiiiii and LSCO^ one 
finds the so-called half- filled stripes, with the density 
of one doped hole per two atoms along the domain 
wall (DW). In contrast, it is clear from a variety of 
experiments that magnetic states within doped Ni02 
planes of the nickelates are filled stripes with density of 



one doped hole per one atom in a DWiSiiiSiLL^ 

The question of filling is not the only difference be- 
tween the nickelate and cuprate stripes, however. Neu- 
tron diffraction measurements performed on Nd-LSCO 
revealed that magnetic peaks are displaced from the an- 
tiferromagnetic (AF) maximum at Qaf — '"'(Ij 1) to the 
points Qj. = 7r(l ± 2e, 1) and Q^, = 7r(l, 1 ± 2e) and the 
shift e depends linearly on hole doping x for x < 1/8, 
while it is almost constant at higher doping. These val- 
ues correspond to a superposition of vertical (01) and 



horizontal (10) DWs. The essentially identical modula- 
tion and doping dependence of e was observed in super- 
conducting crystals of LSCO with x > 0.05. Conversely, 
experiments on LSNO established that spin order is char- 
acterized by the wave vectors Q^ = 7r(l ± e, 1 ± e) with 
e ~ x for X < 1/3, corresponding to a constant charge 
of one hole/Ni ion along a diagonal DW, in agreement 
with the predictions made in the pioneering works by Za- 
anen and Gunnarsson-^ and othersr^iSiiSSiS and empha- 
sized recently,^'* which predicted theoretically the stripe 
order. More precisely, neutron scattering measurements 
have revealed static stripe order in the LNO samples with 
S - 0.105, 0.125, as well as 0.133,^5,26,27,28,29,30,3i ^nd 

even over a wider hole doping regime 0.135 < x < 0.5 in 
the case of LSNOi2i2i^i54S4L& Moreover, the incommen- 
surate (IC) stripe order persists up to x = 0.7 in the 
Nd2-2;Sr2,Ni04 system, ^^ in which the low temperature 
orthorhombic phase seems to extend to a higher dop- 
ing region x < 0.45 as compared to the La compounds, 
where the high temperature tetragonal phase is stabilized 
instead already at a; ~ 0.22>i 

Indications of a charge order (CO) in doped LSNO 
were also found in electron^^ and x-ray diffraction 
studiesi^^iSS*^ Quite recently, the one-dimensional na- 
ture of the stripe modulation within Ni02 planes has 
been directly confirmed in the transmission electron mi- 
croscopy (TEM) studies of charge stripes in LSNOj^I In 
addition, careful examination of the TEM images has 
shown that at a low temperature stripes are mainly cen- 
tered on rows of Ni atoms. However, a mixture of the so- 



called site centered (SC) and bond centered (BC) stripes 
was also observed in some small regions of the sample. 
In contrast to LSCO, charge and spin order in LSNO 
is characterized by the wave vectors Qc = 27r(e, e) and 
Qs — 7r(l ± e, 1 ± e), see Ref. i38j, respectively, with e ~ x 
corresponding to a constant charge density of one hole/Ni 
ion along the diagonal stripe. Note that the doping de- 
pendence of e is exactly the same as that found in the 
seminal Hartree-Fock (HF) calculations within the Hub- 
bard model,i9,20i2i,22,23 

Experimentally the data give a clear evidence that the 
stripe order is most stable at hole doping x — 1/3. In- 
deed, T^Q and T^ increase linearly with x, reach a max- 
imum at a; = 1/3 with 240 and 180 K, respectively, and 
then decrease monotonously upon further doping. The 
particularly robust stability of stripes at x = e — 1/3 
stems from the cooperative spin and charge modula- 
tion expressed by the coincidence of the spin superlattice 
peaks at the wave vectors (l/2±e/2, l/2±e/2) with those 
of CO, given by (e, e). One should also note that e starts 
to deviate gradually from the value given by the e ~ a; 
law above x — 1/3. In fact, CO is even more stable than 
the stripe order at x = 1/2 and forms a checkerboard 
pattern below the transition temperature Tqq ~ 480 K. 
Remarkably, with decreasing temperature, a stripe CO 
sets in below T^q ~ 180 K and its incommensurability 
is twice as large as that of the spin order with the much 
lower onset temperature T^ ~ 80 K^^ The low temper- 
ature competition of the checkerboard and stripe order 
at X = 1/2 has also been clearly indicated by the mea- 
surements of Raman and optical conductivityi^SiiMi In- 
terestingly, above this doping, incommensurability tends 
to saturate rapidly with the value e ~ 0.44.'^^ 

Apparently, CO itself induces commensurate values of 
e and such a commensurability effect seems to be an in- 
trinsic property of the stripe order. At the same time, 
the increased hole density at the DWs suppresses the 
superexchange energy gain, best optimized when all the 
holes are accommodated within stripes. Therefore, it is 
the AF order that drives e to approach the value given 
by the linear relation e ~ x below Tn. Obviously, when 
X — 1/3, both effects cooperate which results in the 
locked-in value of e ~ 1/3 over the entire T < T^q range. 

Furthermore, the temperature dependence of the spe- 
cific heat Cv shows a distinct anomaly at the same tem- 
perature Tqq = 240 K suggesting a CO transition. ^^ This 
conjecture is supported by the temperature dependence 
of the logarithmic resistivity. Indeed, as expected for this 
transition, logp of the sample with x = 0.3 exhibits a 
steep increase precisely at 240 KM. Charge fluctuations 
around T^ = 240 K also lead to conspicuous changes 
in the optical conductivity spectra a{Lo)M- Namely, at 
a high temperature T > Tqq, only a broad peak is ob- 
served with a finite low-energy spectral weight. However, 
when T is decreased down to T^q j the low-energy weight 
below 0.4 eV is gradually transferred to higher energy so 
that the opening of the charge gap is clearly observed. 

A special character of a; = 1/3 as well as a; = 1/2 



doping, is best seen in the magnetic susceptibility x and 
logarithmic resistivity logp, recorded at 200 K, showing 
distinct anomalies at these doping levels4i In this con- 
text it is important to discuss a peculiar behavior of the 
Hall coefficient Rh and the Seebeck coefficient S at two 
representative values of temperature below and above 
Tl% = 240 K, i.e., at 300 K and 210 K. The Hall co- 
efficient at 300 K is almost independent of doping and 
takes a small positive value corresponding to the order 
of one hole carrier per Ni sitei^ The Seebeck coefficient 
shows a similar nearly constant behavior taking, however, 
a negative value. In contrast, below T^q, both Rh and 
S change their signs from negative to positive at a; = 1/3 
and their absolute values are larger than those at 300 
K. In addition, for samples with x = 0.3 and 0.33, Rh 
keeps decreasing with decreasing T, so that the number 
of carriers per Ni site gets reduced even down to 0.01. 

These results indicate that the deviation of x from 1/3 
can be considered as an electronlike {x < 1/3) or holelike 
(a: > 1/3) carrier doping into the e — 1/3 charge-ordered 
insulator with three Ni sites in the unit cell. Hence, for 
the doping level x = 1/3 there is exactly one hole per unit 
cell and such a state is robust and may be considered as 
a half-filled one. Moreover, it would certainly retain this 
feature if the incommensurability had followed precisely 
the relation e = x. However, e has a tendency to shift 
towards 1/3, for both sides of the x = 1/3 point, which 
has important implications for the sign of Rh- On the 
one hand, when x is less than 1/3, the number of holes is 
insufficient for filling up the mid-gap states entirely, i.e., 
the states inside the charge transfer (CT) gap induced 
by stripe order, and in this case the mid-gap states con- 
tain some electrons which become carriers. On the other 
hand, for x larger than 1/3, the number of electrons is 
insufficient for filling up the lower Hubbard band entirely, 
which initially contains holes. Consequently, Rh is ex- 
pected to have the opposite sign to that in the case of 
X < 1/3 (hole carriers). 

Next, the chemical potential shift A/z in LSNO for 
X < 1/3 deduced either from x-ray (XPS) or ultravi- 
olet photoemission spectroscopy^ is suppressed. Cer- 
tainly, this phenomenon cannot be explained within a 
simple rigid-band framework in which /i is expected to 
shift downwards with increasing hole doping. In fact, an 
increase of x in a system with a spatially uniform hole dis- 
tribution should enhance the average hole-hole repulsion 
which, in turn, would result in a higher energy required 
to add one hole to the system, i.e., in a larger \fj,\. There- 
fore, the absence of A/i implies that the average hole-hole 
repulsion remains nearly unaltered upon doping. Such a 
behavior might be easily explained within a stripe picture 
in which a constant hole density at the DWs implies that 
the interwall distance decreases and e increases linearly 
with increasing x. Moreover, a similar suppression of the 
shift has been found below x < 1/8 in LSCO, suggesting 
that an inhomogeneous charge distribution is a common 
feature of both systems. 

Several methods have been employed to investigate the 



stripe phases which go beyond the HF approximation, 
such as: density matrix renormahzation groupp^S slave- 
boson approximation(^2iSSi^ variational local ansatz 
approximationfSS exact diagonalization (ED) of finite 
clusters)^ analytical approach based on variational trial 
wave functions within the string picturej^^^ dynamical 
mean field theory;^ cluster perturbation theory^S and 
quantum Monte Carlo simulations i^SiS They all address 
the crucial role of a proper treatment of local electron cor- 
relations in stabilization of the half-filled stripe phases in 
the cuprates. In spite of this huge effort, it remains un- 
clear whether DWs are centered on rows of metal atoms 
as in SC stripes, or if they are centered on rows of oxygen 
atoms bridging two metal sites in BC stripes. 

Although no evidence was presented yet, it seems that 
the degeneracy of 3d orbitals plays an important role in 
stabilizing filled stripes in the nickelates. In the sim- 
plest picture developed for the cuprates, the Cu'^"'" ions 
forming DWs are spinless, while the Cu^^ of the AF do- 
mains carry a spin S = 1/2. In the nickelates, filled 
DWs are formed of the Ni^+ ions {S = 1/2), whereas 
the AF domains consist of Ni^+ ions (5* = 1). There- 
fore, a realistic Hamiltonian for LSNO has to contain, 
besides the x'^ — y^ orbital which is occupied by one hole 
in the parent compounds of the superconducting cuprates 
(such as La2Cu04 and YBa2Cu306), as well the 3z^ — r^ 
orbital at each ion, so as to account for the high spin 
{S = 1) state of the stoichiometric compound. Indeed, 
when extending their approach to a more realistic four- 
band Peierls-Hubbard model for Ni02 planes, Zaanen 
and Littlewood have shown^ that doped holes prone to 
form diagonal DWs, centered on rows of Ni atoms, with 
a tendency to have a ferromagnetic (FM) alignment of 
the reduced Ni magnetic moments at a DW. In fact, sub- 
sequent multiband HF calculations emphasized the rele- 
vance of the electron-lattice coupling — depending on its 
strength one can obtain either metal- or oxygen-centered 
struct ures i^^i^^ Therefore, it appears that the semiclas- 
sical theory captures the essence of stripe physics in the 
nickelates. 

The stripes were also obtained in the theory using ei- 
ther the HF approximation applied to a multiband CT 
model, in which both the nickel and oxygen degrees of 
freedom are explicitely taken into account jffiiSiiS or by 
the ED of finite clusters within the effective two-band 
model (E3) extended by the coupling of eg electrons to 
the latticeiS However, due to a large number of basis 
states, the latter calculation has solely been done for an 
eight site cluster allowing only for investigating stripe 
phases with a small unit cell observed experimentally at 
high doping levels x = 1/3 and x — 1/2. In contrast, 
the HF approximation used here allows one to describe 
charge and spin modulation with larger unit cells and 
hence it should provide an answer to the important ques- 
tion whether the description of Ni02 planes by the sim- 
plified effective eg model described below yields results 
consistent with the predictions of the CT modelj^2*SiiSS 
An important aspect of these studies, qualitatively dif- 



ferent from the pioneering work on the cuprates where 
the relevant models yield DWs with nonmagnetic Cu'^"'" 
Jgjjg^i9j20j2ij22j23 jg ^]-^g^|. jy^^ j^ ^j^q nickelatcs are formed 

of Ni'^+ ions carrying a finite spin 5=1/2 leading to a 
FM polarization around them. 

In this work we investigate various stripe phases result- 
ing from an effective Hamiltonian for eg electrons relevant 
to the Ni02 layers, and compare them with those ob- 
tained in the doubly degenerate Hubbard (DDH) model. 
We use the HF approximation in order to gain the qual- 
itative insight into possible stripe phases and their sta- 
bility. This approach is the first step to make in order 
to identify possible generic instabilities towards stripe 
phases in Ni02 planes, in analogy to the early work on 
Cu02 planes which predicted the existence of stripes in 
the high Tc cuprates.^® The HF method is well suited 
to compare the energies of different types of magnet- 
ically ordered phases, as it approaches the same limit 
at [/ ^ oo as better mean- field theories, such as slave- 
boson approach or the Gutzwiller ansatzfM In addition, 
the HF approach allows one for efficient calculations on 
large clusters at low temperatures which are necessary to 
obtain unbiased results concerning the stability of vari- 
ous stripe phases, particularly when the magnetic (and 
charge) unit cells are large at low doping. 

The paper is organized as follows. The models with 
orbital degeneracy are introduced in Sec. ^ The energy 
of various stripe phases, filled vs. half-filled, diagonal vs. 
vertical, together with the influence of the various pa- 
rameters of the model, are presented in Sec. IIIII Their 
structures are further described in Sec. IIVI together with 
the mechanism which leads to their stability, illustrated 
by the double occupancy distribution and by the respec- 
tive densities of states. The paper is summarized in Sec. 
Ivl where we also present our main conclusions. 



II. TWO-BAND HUBBARD MODELS 

Even though doped nickelate LSNO is isostructural to 
its cuprate counterpart LSCO, the electronic degrees of 
freedom in LSNO are more involved. In fact, a minimal 
realistic Hamiltonian for LSNO must contain, besides the 
x^ — y^ orbital states included in the cuprate oxide mod- 
els, also the ?>z'^ — r^ orbital states. Such a model of 
interacting eg electrons in a 2D (a, b) plane may be writ- 
ten as follows, 



H — iJkii 



H, 



int 



H, 



(1) 



with two orbital flavors, \x) = x^ — y^ and \z) = 3z^ — r^, 
forming a basis in the orbital space. The kinetic energy 
is described by 



(ij) alia 



with 
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where t stands for an effective (dda) hopping matrix el- 
ement due to the hybridization with oxygen orbitals on 
Ni— O— Ni bonds, and the off-diagonal hopping tf^ along 
a and b axis depends on the phase of the \x) orbital along 
the considered cubic direction. The electron-electron in- 
teractions contain only on-site terms, which we write in 
the following form. 

Hint = ^2^ ("ixT^-ixi +"-izT"'Jzi) 
i 

+ {U -^^Jh)Y^ n^x^iz - 2 Jff Yl ^« ■ ^i^ 

i i 

i 

where U and Jh stand for the intraorbital Coulomb and 
Hund's exchange elements. We also used n.ia — J2a ''^iao- 
for total electron density operators, given by the sum of 
densities in orbitals a — x, z. The last term Hz describes 
the uniform crystal-field splitting between \x) and \z) or- 
bitals along the c axis. 



Hz^lEzYin 






(5) 



The splitting between the eg orbitals originates from 
the tetragonal Jahn- Teller distortion of the NiOg octa- 
hedron. In La2Ni04, however, the octahedron where the 
Ni— O— Ni in-plane (out-of-plane) bond lengths are 1.95 
(2.26) Aj^ respectively, is much less distorted as com- 
pared to 1.89 and 2.43 A bonds in La2Cu04r^ which 
reflects the difference in electron filling. In what follows 
we consider only a realistic positive Ez favoring, due to 
elongated octahedra, the \z) electron occupancy over the 
I a;) occupancy in doped compounds. 

In order to elucidate the mechanism leading to the 
stabilization of the stripes, we also investigate the DDH 
model. It is obtained using the hopping matrix. 



tap 



1 
1 



(6) 



instead of Eq. (PJ. Note, that the total bandwidth is 
W = Qt in the Cg model and M^ = 4t in the DDH model 
(at i?2 = 0), but the average diagonal hopping elements 
are the same — this ensures that electrons are approxi- 
mately similarly correlated for a given U in both cases, 
while the electron-electron interactions have in the DDH 
model the same form given by Eq. Q), provided one la- 
bels the orbitals as x and z. 

The model given by Eq. Q is very rich. Even when 
restricting oneself to solutions with only two atoms in the 
unit cell, one obtains an intricate competition between 
numerous FM and AF phases, even on the level of mean 
field (MF) approximation.^^ In particular, we obtained 
phase diagrams showing, in the regime away from half 
filling, a clear tendency towards ferromagnetism for Jh ~ 
C//4, or at least to the C-AF phase for Jr — 0.15C/. 
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FIG. 1; (Color online) Fermi surfaces of the higher (top) and 
lower (bottom) bands, as obtained for different values of the 
chemical potential fi in the two-band tight-binding model for 
degenerate Cg orbitals {E^ — 0). Note a different scale for 
both Fermi surfaces. 



Some of the richness of the present model (Q) can be 
traced back to the peculiar features of the electronic dis- 
persion for Bg electrons in two dimensions^ Unlike for 
spins, two bands obtained by diagonalizing the Hamil- 
tonian (^) in the noninteracting limit are nondegenerate 
(except for some high symmetry directions)!^ 



£k,. = -t{X + Y) 



vty/X^ + Y^-XY- e(X + Y)+ e^, (7) 



with 1/ = ±1, X = cos/ct, Y 



cos ky, and 



El 

2t 



(8) 



For e = the corresponding Fermi surfaces are depicted 
in Fig. n At half filling one piece of the Fermi surface is 
building a " Swiss cross" , while the second one is building 
a square. The filling of the upper band is then by 1/4 
electron, while that of the lower band is by 3/4 electron 
per spin. The Fermi surface of the upper band shrinks 
rather fast with increasing doping x and below fj, = —t 
{x = 0.93) one finds that this band is empty and all 
electrons are in the lower band. However, both bands 
are partly filled in the doping regime x < 0.4 relevant for 



the nickelates, and they will both play a role also when 
electron interactions are included. 

At half filling the edges of the pieces of the Fermi sur- 
face (FS) for both bands are connected by the nesting 
vectors: Qi — (0, tt) and Q2 = (tt, 0). At the same time, 
the edges of one piece are connected to the edges of the 
other one by the nesting vector Q3 = (tt, tt). In addition, 
at half filling both components of the velocities vanish on 
their respective Fermi surfaces. Therefore, several com- 
peting instabilities are expected once the interaction is 
turned on. Their relative importance can be estimated 
by computing the spin-spin correlation function in the 
noninteracting limit for the corresponding wave-vectors. 
In the static limit it reads: 



xs(Q) 



2 >r^ y^ /F(ek^) 



fp i^k+Q^) 



N 



£k+Qiy — £k;i 



[/„^(k)[/t^(k)(7^,(k -f Q)[/tjk + Q) .(9) 



Here the matrices C/(k) diagonalize the hopping matrix 
Eq. (O, and £k^ are given by Eq. {71). 

Since ferromagnetism is an expected feature of the 
model (Q), we also evaluate the homogeneous spin-spin 
correlation functions. It shows an intricate interplay of 
all the three instabilities (see Fig. EJ. They were all re- 
ported, regardless of the approach, see for instance Refs. 
163 and|63 One fi nds th at the van Hove singularity takes 
place at /x = — vo + e^t for the crystal-field splitting e, 
see Eq. ^. At this point the line where the velocity 
vanishes has a large component on the FS, and the sin- 
gularity is logarithmic. Note that there is no nesting 
vector connecting different pieces of the FS, as can be 
understood by looking at Fig. ^ Here one finds a dis- 
tinct difference between the present Cg model and the 
DDH model, where instead the FS is nested and takes 
the shape of a square for both bands at the van Hove 
point (for a; = 0). Therefore, nesting and van Hove sin- 
gularities in the eg electron model influence the physical 
properties at different electron densities than in the DDH 
model. 




FIG. 2: Static spin-spin correlation functions Xs(Q) for 69 
electrons for increasing doping, as obtained for: (a) E^ — 0, 
and (b) E^ — t. Different lines correspond to: Q — (0, 0) 
(solid line), Q = (0,7r) (dotted line), and Q — (tTjTt) (dashed 
line). Van Hove singularity for Q — (0,0) is obtained at 
doping higher than that realized experimentally in nickelate 
LSNO oxides: (a) x ~ 1.49, and (b) x ~ 1.68. 



In the following the Hamiltonian given by Eq. Q is 
treated in the HF approximation. We introduce explicity 
four "diagonal" expectation values, (riiaa), and four "off- 
diagonal" ones, (cJ^o-Cja^), with a, P = x,z and a ^ 
/3, that we determine self-consistently in order to get a 
complete characterization of the ground state. Due to 
symmetry the last four Fock averages vanish for a 2D 
square lattice, and one gets the simpler MF Hamiltonian: 



HmF - 22 '^iac7^V,al3,<TCj/3„, 
ij.af3.(7 



with. 
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(11) 



M,,„ = i J2 [i^U - 5JH){n,) -XAU + JH){m,) 

i 

- ([/ - 5J„){o,) - \,{U - JH){P^)] + \E,, 
M^za = 3 E [(3f^ - 5Jff)(n,;) - X,{U + JH){m,) 

i 

+ (U- 5Jh){o^) + \a{U - JH){p^d] - \E,, (12) 

where we introduced A^ = ±1 for a =] (J,) spin and 
Aq = ±1 for a = x(z) orbital, and the operators: 



aa 

aa 

aa 
Pi = / ^ XqAgniaa, (13) 



ra 



for total charge and magnetization density, and for or- 
bital anisotropy in the charge and magnetization distri- 
bution, respectively. 



III. GROUND STATE ENERGIES 

Site centered and bond centered stripes represent 
equally relevant candidates for the ground state. How- 
ever, it numerically turns out that the energy of SC 
stripes is substantially larger than the one of the EC 
stripes, in agreement with the early study of Yi et dlm^ 
(see Table UTTl in Sec. I^as well). We therefore concen- 
trate on BC stripes. They are made out of pairs of atonrs 
with a FM spin polarization within the model Hamil- 
tonian (^ in the wide doping regime 0.05 < x < 0.4. 
In order to obtain unbiased results, we performed calcu- 
lations on large clusters implementing the symmetry of 
stripe phases in reciprocal space. We worked on squared 



clusters with the hnear dimension along the x direction 
chosen as an even multiplicity of the elementary stripe 
unit cell dimension. Their sizes are listed in Table HI 

To ensure that the model (^ is indeed relevant for the 
nickelates, it is necessary to focus on appropriate values 
of parameters U, Jh, and E^. First of all, x-ray absorp- 
tion (XAS) measurcments^S as well as from XPS and 
bremsstrahlung-isochromat spectroscopy studies^i of the 
electronic structure of LSNO suggest that LSNO is a CT 
insulator with nearly the same CT energy A as that of 
NiO. Therefore, the same parameters could be accepted 
as used before in the self-consistent Born calculationsji^ 
which reproduced quite well photoemission spectra of 
NiO. Secondly, since our approximation has a tendency 
to overestimate the effect of the Coulomb repulsion in 
LSNO, we used a somewhat smaller value of the Coulomb 
repulsion U than that adequate for NiO, i.e., we set 
[/ = A = 5 eV. Next, as the in-plane Ni-O-Ni bond 
length in La2Ni04 of 1.95 A(see Ref. |65|) is very much 
the same as the shorter bond of 1.89 A in La2Cu04fS& 
we set the hopping amplitude tpd between the pa or- 
bitals and \x) ~ \x'^ — y^) orbitals to be as in LSCO, 
i.e., tpd — 1.47 eVi2^ This in turn yields an effective in- 
plane Ni-Ni hopping i^^ = {t^^/A) = 0.43 eV. However, 
it is more convenient to take the effective (dda) hopping 



0.00 



element connecting two \z) 



3z^ 



orbitals along the 



c-axis as the energy unit t, related to t^^ via the Slater- 
Koster relation t = 4i^^/3 ~ 0.6 eV, so that U ~ 8t. 

The value of Hund's exchange between t2g electrons in 
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FIG. 3: Free energy gains SF ^ of the VBC (top) and DBC 
(bottom) stripe phases with respect to the AF phase H15|l as 
functions of doping x, obtained at temperature pt — 100 for 
the Eg model. Cluster sizes, distances between DWs d, and 
parameters are the same as in Tabled] 



NiO was estimated as J^ 



0.8 evil It is related to Jh ^°' '^ '^''*™^' ^^^^^^ '^' '™P^' ''^^*^°^ 

Jh = Jh + ^i 



75 



(14) 



TABLE I: Comparison of the ground state free energy F for 
the VBC (left) and DBC (right) stripe phases as found in 
the Eg model for increasing doping x. Stripes are separated 
by d = 3, . . . , 11 lattice constants in clusters of size L x L. 
Parameters: U = 8t, Jh = 1.5t and E^ — 0. 







VBC 






DBC 




X 


d 


Lx L 


F/t 


d 


Lx L 


F/t 


0.05 


11 


88x88 


2.8811 


11 


88x88 


2.8794 


0.06 


11 


88x88 


2.8376 


11 


88x88 


2.8326 


0.07 


11 


88x88 


2.7946 


11 


88x88 


2.7870 


0.08 


11 


88x88 


2.7519 


11 


88x88 


2.7430 


0.09 


10 


80x80 


2.7099 


10 


80x80 


2.7001 


0.10 


10 


80x80 


2.6677 


9 


72x72 


2.6570 


0.11 


9 


72x72 


2.6258 


8 


64x64 


2.6140 


0.12 


8 


64x64 


2.5839 


7 


84x84 


2.5712 


0.14 


7 


84x84 


2.5001 


6 


72x72 


2.4853 


0.16 


6 


72x72 


2.4162 


6 


72x72 


2.3997 


0.18 


5 


80x80 


2.3329 


5 


80x80 


2.3130 


0.20 


5 


80x80 


2.2485 


4 


80x80 


2.2291 


0.25 


4 


80x80 


2.0394 


4 


80x80 


2.0155 


0.30 


3 


72x72 


1.8325 


3 


72x72 


1.8029 


0.40 


3 


72x72 


1.4619 


3 


72x72 


1.4081 



where B stands for a Racah parameter^ Taking into ac- 
count that B ~ 0.13 eV for Ni02r^ one finds Jh = 0.93 
eV, i.e., Jh — 1.5t. Indeed, it has been shown that Jh = 
1 eV reproduces the experimental band gap and the mag- 
netic moment of La2Ni04.^^ Finally, we set E^ ~ t — 0.6 
eV as a realistic value of the crystal field splitting in 
the nickelates. On the one hand, band structure calcula- 
tions in the local density approximation predict the crys- 
tal field splitting between eg orbitals to be 0.5 eVm^ On 
the other hand, XAS spectra reveal a somewhat larger 
splitting of 0.7 eV,*- a value also deduced from the opti- 
cal spectroscopy,"^^ 

The performed calculations demonstrate a robust ten- 
dency towards stripe formation expressed by negative 
free energy gains of the vertical bond-centered (VBC) 
and diagonal bond-centered (DBC) stripe phases with 
respect to the AF phase. 



6F = Fs~ Fai 



(15) 



for the Cg model 



(E3). This energy gain \SF\ increases 
with doping x, as shown in Fig. 13 The 6F curves cross 
each other for decreasing distances between DWs d which 
demonstrates the tendency to the gradual formation of 
stripe phases with smaller unit cells upon increasing dop- 
ing. Thereby, diagonal structures are significantly lower 
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FIG. 4: Free energy gains 5F as in Fig. |3 but with finite 
crystal field splitting Ez = t. Cluster sizes, distances between 
DWs d, and parameters are tlie same as in Table |I] 
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FIG. 5: Free energy gains 5F of the VBC (top) and DBC 
(bottom) stripe phases with unit cells of length d, as a func- 
tion of doping X, obtained for the DDH model. Parameters: 
U = 8i, Jh = 1.5f, E^ = 0, and (31 = 100. 



in energy than vertical ones for a given fixed doping, but 
especially in the large doping regime x ~ 0.4. The ro- 
bust stability of the DBC stripe phases with respect to 
the VBC ones is illustrated more transparently in Ta- 
ble m where we compare the ground state free energy 
F for both structures. Note that a similar variation of 
the optimal distance d between DWs suggests the same 
optimal stripe filling. 

The main effect of the crystal field splitting is to reduce 
the stability of both stripe phases with respect to the uni- 
form AF phase. This appears to be rather puzzling at 
first sight. One finds, however, that realistic positive E^ 
promotes the \z) orbitals with a narrow band. On the one 
hand, the electron distribution induced by finite E^ sup- 
presses substantially the AF superexchange energy gain 
which becomes J^z = '^'tlz/U rather than J^x — 'it'^x/U- 
On the other hand, one would expect that such charge 
redistribution strongly increases the kinetic energy gain 
due to a wider band accessible for holes, propagating es- 
pecially easily along DWs where the AF order is partly 
suppressed. Nevertheless, this energy gain is easily over- 
compensated by the kinetic energy loss due to hopping 
perpendicular to the stripes. In addition, our studies 
of the BC stripe phases within the single-band Hubbard 
model have shown that the largest kinetic energy gain 
is released on the bonds connecting pairs of ferromag- 
netically coupled atoms located in the DWs.^'' Note that 
the FM order of the DW spins is substantially stabilized 
by the off-diagonal txz hopping in the Cg model, yield- 
ing low-energy charge excitations that lead to the FM 
superexchange, Jxz = 'it'^z/iU — 3Jh)- Altogether, when 



one orbital is sufficiently favored by finite crystal field 
over the other one, these low-energy processes are effec- 
tively blocked, explaining the enhanced stability of the 
AF order with respect to the BC solutions (c/. Figs. |2| 
and ID). 

A further qualitative point concerns the influence of a 
finite crystal field splitting Ez between the |a;) and \z) 
orbitals on the stability of DW structures. As depicted 
in Fig. ^ a more realistic value E^ = t seems not to 
promote noticeably any stripe phase over another one 
and one still recovers DBC stripe phases as the ground 
state. We therefore conclude that it is not the crystal 
field Ez that is responsible for a different orientation of 
DWs in the nickelates from that observed in the cuprates. 
Later on we shall see that finite Ez has also only a little 
visible effect on optimal stripe filling. 

Let us now shortly return to the experimental char- 
acterization of the stripe phases in LSNO. Their peri- 
odicity at a given doping x is described by both charge 
Qf, = ±7r(2e, 2e) and spin Q^ = ■7r(l ± e, 1 ± e) wave 
vectors along the direction diagonal with respect to the 
Ni— O bond. The incommensurability e corresponds to 
the inverse of the distance d between DWsfii2i2i^iS*SiL& 



(16) 



Further, e first increases continuously with doping x and 
follows the linear relation e = a; in the wide doping range 
X < 1/3, while next it gradually saturates to the value e ~ 



8 



0.44.^ Such a behavior of e indicates a fixed hole density 
of one hole per Ni ion in a DW, and is consistent with 
the HF predictions both in the single and in multiband 
niodels.^^i^°i^^i^^i^'^i^°i"i^^ Finally, a fixed hole density 
along the DWs results in the pinning of the chemical 
potential ^ for x < 1/3, whereas a large (~ 1 eV/hole) 
downward shift appears in the higher doping reginie4& 

It is interesting to establish whether the above results 
concerning the stripe phases in the ground state and their 
variation under increasing doping appear solely in the 
realistic e^ model, or they are generic and remain also 
a common feature of the DDH model with only diag- 
onal hopping elements given by Eq. I^. To allow a 
meaningful comparison, we set the same values of micro- 
scopic parameters as we those chosen for the eg model, 
i.e., U = 8i, Jh = 1.5t, and we consider only the E^ = 
case imposed by symmetry between the equivalent or- 
bitals. It restricts the solutions of the DDH model to 
(o) = (p) = (see below). 

One finds that the free energy gains SF of the VBC 
and DBG stripe phases with respect to the AF phase for 
the DDH model (Fig. O are qualitatively similar to those 
of the eg band model. First of all, for a fixed doping x 
diagonal stripe structures are again significantly lower in 
energy than vertical ones (c/. also Table lll| . Secondly, 
also in this case, we recover a gradual crossover towards 
stripe phases with smaller unit cells upon increasing dop- 
ing. Note, however, that in contrast to the predictions 
made within the Cg model, structures with vertical (diag- 



TABLE II: Comparison of the ground state free energy F 
per site for the VBC (left) and DBC (right) stripe phases as 
found in the DDH model for increasing doping x. Stripes 
are separated by d lattice constants in clusters of size L x L. 
Parameters: U = 8i, Jh ~ l-5t, and Ez = 0. 







VBC 






DBC 




X 


d 


LxL 


F/t 


d 


LxL 


F/t 


0.050 


11 


88x88 


3.1095 


11 


88x88 


3.1000 


0.055 


10 


80x80 


3.0912 


11 


88x88 


3.0804 


0.060 


9 


72x72 


3.0729 


11 


88x88 


3.0610 


0.065 


8 


64x64 


3.0547 


10 


80x80 


3.0417 


0.070 


8 


64x64 


3.0364 


9 


72x72 


3.0225 


0.080 


8 


64x64 


3.0001 


8 


64x64 


2.9839 


0.090 


7 


84x84 


2.9636 


7 


84x84 


2.9454 


0.100 


6 


72x72 


2.9269 


7 


84x84 


2.9070 


0.110 


5 


80x80 


2.8902 


6 


72x72 


2.8683 


0.120 


4 


80x80 


2.8538 


5 


80x80 


2.8301 


0.140 


4 


80x80 


2.7807 


5 


80x80 


2.7530 


0.160 


3 


72x72 


2.7077 


4 


80x80 


2.6760 


0.180 


3 


72x72 


2.6347 


4 


80x80 


2.5996 


0.200 


3 


72x72 


2.5626 


3 


72x72 


2.5227 


0.300 


3 


72x72 


2.2129 


3 


72x72 


2.1531 


0.400 


3 


72x72 


1.8809 


3 


72x72 


1.8116 



onal) DWs separated by a distance d > 5 are the lowest 
energy solutions only in a narrow doping regime x < 0.12 
(x < 0.15), respectively. 

For a complete characterization of stripe phases we 
define 



Nh 



LNuw' 



(17) 



where N^ is the number of doped holes and iVow is the 
number of domain walls. Our findings concerning the 
properties of BC stripe phases are summarized in Fig. 
showing the doping dependence of the incommensurabil- 
ity e, the stripe filling v, and the chemical potential /i 
for both the VBC (left) and DBC (right) stripe phases. 
These quantities for the respective ground states are de- 
duced from Figs. |2| El ^^^d |S1 In each case the data 
points correspond to the middle of the stability region of 
the stripe phase of a given type. The only exception are 
the d = 3 phases — for them e, v, and /i, are plotted for 
the actual minimum of the free energy F. 

In agreement with the experimental data for LSNO, 
indicated here by a solid line in panels (a) and (b) of Fig. 
El one observes that e follows the law e~a;uptoa;~0.2 
and then it tends to saturate to the highest possible value 
for the BC stripe phase, i.e., e — 1/3, as there are no 
BC stripe phases in which DWs are separated by d = 2 
lattice spacings. Further, in the regime where e follows 
linearly the doping level x, decreasing stripe periodicity 
allows the system to maintain nearly fixed filling i/ [cf. 
Figs, (njc) andEJd)], pinning simultaneously the chemical 
potential /i, as shown in Figs, ^c) andlH^f). Remarkably, 
the ground state of both the VBC and DBC stripe phases 
is characterized by the optimal filling 0.9 hole/Ni, very 
close indeed to the experimental value one hole/Ni ion, 
and the optimal filling remains almost unaltered in the 
model with a finite crystal field E^ = t. 

Regarding the chemical potential shift A/i in the dop- 
ing regime x > 0.2, one finds that it exceeds the exper- 
imental value ~ —1.0 eV/hole nearly by a factor of 2. 
Indeed, assuming the effective hopping t = 0.6 eV, one 
obtains A/i ~ -2.2 (-2.1) eV/hole for S^ = {E^ = t), 
respectively. Therefore, we conclude that the present 
effective model can only explain qualitative trends and 
one needs to carry out calculations within more realistic 
multiband models with oxygen orbitals included explic- 
itly in order to obtain quantitatively the experimental 
data. 

Let us verify now whether the established results con- 
cerning the doping dependence of e and v are also ob- 
tained within the DDH model. As the structures with 
vertical (diagonal) DWs separated by a distance d > 5 
appear in the DDH model only in a narrow doping regime 
X < 0.12 {x < 0.15), respectively, one observes here a fast 
variation of the optimal distance d which should result in 
a small optimal stripe filling. Indeed, as depicted in Figs. 
|SJa) andEJb), the low-doping part of e, being linear in x, 
has a larger slope which exceeds the experimental value 
in LSNO roughly by a factor of 2 (1.5) in the case of the 
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FIG. 6: (Color online) Doping dependence of stripe phases for 
the VBC (left) and DBC (right) ground states: (a) and (b) 
magnetic incommensurability e; (c) and (d) stripe filling v\ (e) 
and (f) chemical potential [i, (in the units of f), as deduced 
from the data of Figs. |3 0| and |^ Filled (open) symbols for 
eg (DDH) model; circles and squares for Ez = and Ez — 
t. Solid lines in panels (a) and (blshow the experimental 
behavior of e in LSNO given in Ref. |3- 



VBC (DBC) stripe phase, respectively. Consequently, 
the optimal stripe filhng in the former case is substan- 
tially reduced down to j/ ~ 0.55 and in the latter case 
— down to J/ ~ 0.65 [see Figs. EJc) andlHId)]. Finally, 
one finds that the chemical potential starts to decrease 
rapidly at x ~ 0.14 in the VBC stripe phase [Fig. EJe)] 
and at slightly larger doping x ~ 0.16 in the DBC phase 

[Fig. Elf)]. 




FIG. 7: (Color online) Charge and magnetization distribu- 
tion in the filled VBC (left) and DBC (right) stripe phase 
found in either the tg model (filled circles), or in the DDH 
model (open circles) on a 64 x 64 cluster at doping x = 1/8: 
local hole densities nh{lx) (top row); modulated magnetiza- 
tion densities 771-^ {la:) (second row); local hole Oh{lx) (third 
row); and local modulated magnetic Ptt{Ix) orbital polariza- 
tion (bottom row). For more clarity, the data points for the 
DDH model are shifted by one lattice constant from the origin 
of the coordinate system. Parameters: U = 8t, Jh = 1.5i, 
Ez = 0, and pt = 100. 



IV. MECHANISM OF STRIPE FORMATION 



A. Charge and magnetization distributions in the 
stripe phases 

In order to find out the reason of such a vast dis- 
crepancy between the predictions made in the Cg and in 
the DDH model, let us now investigate closer the prop- 
erties of idealized filled and half-filled DWs: half-filled 



vertical bond-centered (HVBC) and half-filled diagonal 
bond-centered (HDBC). A complete characterization of 
the charge and magnetization distributions for a system 
with orbital degeneracy is given by: the local hole density 
x{lx), local modulated magnetization density mTr{lx), lo- 
cal hole orbital polarization Oh{lx), and local modulated 
magnetic orbital polarization Pt^{Ix), defined using the 



10 



operators given in Eqs. (|13|) as follows: 
nfiilx) = y^^nhajlx), 

a 

m-^{lx) = (-1)'" ^ Xanaa{lx), 

cxa 

a 
PttQ-x) = (-1)'"- y^ \a\ana„{lx). 



(18) 



Here naCT(^a;) is the local orbital charge density for spin 
a, whereas rihai^x) denotes the local orbital hole density, 



^ha [''X 






^aa [''X 



(19) 



To better appreciate the differences between the DDH 
and the Cg model, we compare in Fig. [7| the local hole 
nh{lx) and the modulated magnetization densities mT^{lx) 
of the filled VBC and DBC stripe phases found at tem- 
perature fit = 100 in either model on a 64 x 64 cluster at 
doping X — 1/8 for the standard parameters, i.e., U — 8t, 
Jh = l-5i, and Ez — 0. For completeness we also show 
the local hole Oh{lx) and the local modulated magnetic 
PTri^x) orbital polarization. These quantities are finite 
only in the Cg model, while Oh{lx) = Pn{lx) = by sym- 
metry in the DDH model. Since Oh is positive the holes 
are located in \x) orbitals, and the polarization of the 
two sites on the DWs is FM (alternation of signs in p^r 
at the DWs). 

The observed differences, especially pronounced at DW 
atoms, directly follow from the fact that, a different ef- 
fect helping to reduce double occupancy at those sites 
is effective in each model. Namely, in the eg hopping 
model Q one finds large positive Oh{lx) at DWs, which 
means that it is energetically advantageous to optimize 
the kinetic energy of holes by putting them into the \x) 
orbitals where the hopping clement txx is larger and gives 
a wide band. On the one hand, in the HF the only way 
to optimize the on-site energy is to develop a strong spin 
polarization which, in turn, would noticeably reduce the 
kinetic energy gain. On the other hand, such disadvanta- 
geous suppression can be avoided by a strong reduction 
of the electron density. This explains why the hole den- 
sity nh{lx) along DWs in the eg model is larger, as com- 
pared to the corresponding value found in the DDH case. 
Indeed, in the latter case both bands are equivalent, re- 
sulting in Oh{lx) = PnO-x) = 0. Hence this model yields a 
more localized stripe phase with a larger magnetization 
|w7r(^a;)| at DWs than the one obtained in the eg model 
(c/. Fig.EJ. 

For completeness, in Fig. |Hlwe compare the hole den- 
sity, nh{lx), and the spin density, m.^{lx), profiles of the 
half-filled stripe phases. The unit cells are smaller by 
a factor of two and the amplitude of the charge density 
wave is correspondingly reduced. Even though the over- 
all shape of the DWs looks very much the same in both 
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FIG. 8: (Color online) Charge and magnetization distribu- 
tions in the half-fiUed HVBC (left) and HDBC (right) stripe 
phases found in either the Cg (filled circles) model, or in the 
DDH model (open circles) on a 64 x 64 cluster at doping 
X = 1/8. The meaning of the different panels and shift of 
the DDH data as in Fig. Q Parameters: U = 8f, Jh ~ 1.5t 
E^ = 0, and /3i = 100. 



models, a larger magnetization |m7r(Z2;)| at DWs is found 
again within the DDH Hamiltonian 10) . 



B. Double occupancy distribution 



Furthermore, important information about the nature 
of stripe phases is provided by the averaged intraorbital 
double occupancy at site Ix, 



D{lx) = ^na^{lx)nai{lx), 



(20) 
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FIG. 9; (Color online) Double occupancy distribution in the 
filled VBC (left) and DBC (right) stripe phase shown in 
Fig- 13 intraorbital double occupancy D{lx) (top row), and 
two local interorbital D^" (Ix) (middle) and D%'^{lx) (bottom) 
double occupancies. Filled (open) circles denote the results 
found in the eg (DDH) model, respectively. Parameters and 
filled/empty circles as in Fig. Q 



as well as by two local interorbital double occupancies. 



DTz{lx) = '^nx,y{lx)nzs{lx) 

a 



(21) 
(22) 



where a = —a. In a system with isotropic charge distri- 
bution the double occupancies are site independent, while 
in stripe phases they vary with a characteristic periodic- 
ity following from the size of the stripe unit cell. From 
the structure of the Coulomb interaction Eq. Q one ex- 
pects D < D%'^ < I?^J for a fixed hole density. In a stripe 
phase these double occupancies are locally suppressed by 
the magnetic and orbital polarizations in a way that de- 
pends both on the shape of the stripe, and on the form 
of the one-electron Hamiltonian. Let us first discuss the 
average double occupancies (|2U|) - (|22|I of the stripe phases 
depicted in Fig. |7| Remarkably, in the DDH model (open 
circles in Fig.|3) the on-site energy is predominantly opti- 
mized by the reduction of the high-energy configurations 
with opposite spins which lead to intraorbital D(lx) and 
interorbital D^'^{lx) double occupancies, so that the sys- 
tem might create a large number of DWs in the unit cell 
even in the low doping regime, consequently reducing 
the optimal stripe filling. Such a strong reduction of the 



FIG. 10: (Color online) Double occupancy distribution in the 
half-filled stripe phases shown in Fig. |H] HVBC (left) and 
HDBC (right). Different double occupancies as in Fig. |5| 
Parameters and filled/empty circles as in Fig. |H1 



double occupancies for opposite spins follows from the 
diagonal hopping which gives a weaker electronic mixing 
between the sites of opposite spin due to a lower kinetic 
energy ii^ 

In contrast, in the Cg model (filled circles in Fig. O, 
the Coulomb energy is mainly optimized by the reduction 
of the low-energy interorbital D^^{lx) double occupancy 
resulting in a smaller magnetization |?n7r(^2:)| at DWs as 
compared to the one found in the DDH model. How- 
ever, reduced |r7i^(Za;)| at those sites allows the system to 
better optimize the kinetic energy gain which then over- 
compensates a large on-site energy only when the optimal 
filling is close to one hole per Ni site, meaning that for 
a given doping level, the DWs should be separated by a 
larger distance as compared to predictions made in the 
DDH model. On the other hand, the robust stability of 
the DBC stripe phases with respect the VBC ones can 
be understood as following from a stronger reduction of 
all double occupancies (|20|I - H22|I in the former phase. In 
fact, the doubly occupied configurations cannot be then 
excited between the nearest neighbor DW sites, leading 
to a lower total energy. 

In case of half-filled stripes (Fig. I10|l the double occu- 
pancies for opposite spins {D{lx) and D%'^{lx)) are only 
slightly lower than for filled stripes, but the interorbital 
double occupancy for equal spins {D^zi^x)) is signifi- 
cantly higher, reflecting the higher electron filling at DW 
atoms in these phases. Due to the finite Hund's exchange 
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FIG. 11: Densities of states for the DBC stripe phase in the 
Eg model at x = 1/8 doping: (a) partial density of states 
Nx{i^) for \x) orbital, (b) partial density of states N^{lu) for 
\z) orbital, (c) total density of states A'^(a;). Panel (d) shows 
for comparison the total density of states N(uj) as found in 
the DBC stripe phase within in the DDH model. Parameters 
as in Fig. 



considered here, Jh = 0.15t, the high-spin states are pro- 
moted. Altogether one finds similar generic features seen 
above for the filled stripes (Fig.lHl): (i) both double occu- 
pancies for the opposite spins are stronger reduced in the 
DDH model, and (ii) the double occupancy by the same 
spin Z?!JJ ^t the DWs is stronger suppressed in the Cg 
model, resulting from a smaller magnetization at these 
sites than in the DDH model. This second effects seems 
to play a role in stabilizing the HDBC stripe phase in the 
DDH model rather than the DBC phase, found for the 
realistic Cg model. 



termine the DOS projected on the a orbital as: 

^"('^) = ^ E E l*.a.(k)P<5(c. - ek.), (23) 



k i{T 



where '^^^^^{k) are the eigenvectors of the MF Hamil- 
tonian Eq. I|1U|) . It is calculated using a histogram of 
the corresponding eigenvalues. It is apparent from Figs. 
im and El that the spectra consist of several subbands. 
However, two Hubbard subbands are well visible both 
for filled and half-filled diagonal BC stripes: lower Hub- 
bard band (LHB) and upper Hubbard band (UHB). The 
Hubbard subbands arise primarily due to the AF polar- 
ization and the gap between them is proportional to the 
magnetization. 

Panels (a) and (b) in Figs. ITTl and [HI show the DOS 
Na{uj) H23|l . projected on the orbital \x) and |z), re- 
spectively. Here, one finds that the mid-gap bands have 
mainly \x) character close to the gap between them — 
the \x) orbitals optimize the kinetic energy of holes doped 
into them, whereas the vast majority of the more local- 
ized I z) states of DW atoms belongs to the energy regimes 
within the LHB and close to the UHB. 

The total DOS for the DBC stripes in the Cg model, de- 
picted in Fig. IllT c') , reveals mid-gap states which demon- 
strate the one-dimensional character of the transport in 
the stripe phases. Note that due to the AF spin modu- 
lation along the BC DWs themselves, one finds as well a 
distinct gap between two mid-gap bands lying within the 
Mott-Hubbard gap. This gap for the filled DBC stipes 
is lower in the eg model [Fig. Illf c')] than in the DDH 
model [Fig. Illf d)]. as the holes are more dclocalizcd in 
the former case. This holds as well for the HDBC stripes, 
see Figs. I12f c') and ll2f d). We have shown before that the 
magnetization on the DW sites is larger for the HDBC 
stripes than for the DBC ones, and therefore the gaps be- 
tween the mid-gap states are also larger in Fig. 1121 than 
in Fig. ini 

Remarkably, both for the filled and half-filled diagonal 
BC stripe phases of Figs. E] and E| the lower mid-gap 
states are to some extent localized above the LHB in the 
eg model, whereas they are clearly almost merged with 
the LHB in the DDH model. This behavior follows from 
a smaller on- wall magnetization m.^{lx) in the eg model. 
Therefore, in the DDH case only a few holes occupy the 
low-lying mid-gap bands, explaining the reason of sub- 
stantial suppression of the optimal stripe filling in this 
case. Instead, filled stripes are favored in the Cg model 
— they have a larger fraction of holes in the lower mid- 
gap states. Note, however, that also in this case the 
stripes are metallic, with some spectral weight (mainly 
in \x) orbitals) above w = /x. 



Densities of states 



DISCUSSION AND SUMMARY 



Important differences between the two models strongly 
influence the density of states (DOS). Here we first de- 



We summarize our key results concerning the stabil- 
ity of stripe phases in the nickelates. We investigated 
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FIG. 12: Densities of states for the HDBC stripe phase in 
the Eg model a.t x = 1/8 doping, panels (a)-(c) as in Fig. 
111! Panel (d) shows for comparison the total density of states 
N{uj) as found in the HDBC stripe phase within in the DDH 
model. Parameters as in Fig. |H| 



this problem using two Hubbard models — the realis- 
tic model for Cg electrons, and the DDH model which 
frequently serves as the simplest approach to the or- 
bital degeneracyJi It is quite encouraging that the DDH 
model gives already diagonal stripe phases as the most 
stable structures, in contrast to the nondegenerate Hub- 
bard model where the vertical (horizontal) stripes are 
found in the HF approximation. However, one finds half- 
filled BC stripes in the DDH model, and it is only in 
the realistic eg model for doped nickelates that one ob- 
serves a generic tendency to promoting filled stripes over 
the half-filled ones. As an example of this behavior, we 
give the free energies of filled and half-filled stripe phases 
found within both (cg and DDH) models for x = 1/8 
doping in Table Hill First of all, one finds that the ener- 
gies of BC phases are considerably lower for the realistic 
model of eg band due to the higher kinetic energy gains 
which result from the off-diagonal hopping. ^^ Most im- 
portantly, for the parameters relevant to LSNO used for 
the data of Table IIIII one finds that the DBC phase has 



the lowest energy, i.e., indeed experimentally observed 
filled diagonal stripes. In contrast, the SC stripe phases 
with unpolarized DW sites and large double occupancies 
have much higher energies, and the best of them is not 
diagonal, but vertical (horizontal) one. We also note that 
even though the DDH model with two equivalent orbitals 
clearly favors diagonal DWs, it stabilizes instead of filled 
the half-filled diagonal phases. 

As in the cuprates, the coexisting charge and magnetic 
order in diagonal stripe phases in the nickelates is a re- 
sult of the compromise between the kinetic and magnetic 
energies — the magnetic energy is gained in the AF do- 
mains, and the kinetic energy is gained mainly along the 
DWs. The BC stripes are favored as then the magnetic 
energy can be gained f Table llTT|l not only in the AF do- 
mains, but also on the DW magnetic sites. Finally, when 
the DWs are filled, more kinetic energy is gained in the 
eg model while then the off-diagonal hopping is allowed. 
Altogether, this mechanism shows little sensitivity to the 
small crystal field splitting. For moderate values of the 
latter it tends to eliminate stripe phases entirely rather 
than promoting one stripe phase over another one. 

In spite of the remarkable success of the present study 
which gave stable diagonal stripe phases, one has to ad- 
mit that the predicted electronic properties show some 
difference to the experimental ones. In fact, our system- 
atic studies of stripe phases completed within the rel- 
evant model for Cg orbitals, where several phases sep- 
arated by different lattice spacing varying from d — 3 
to d = 11 were considered, have revealed that the opti- 
mal stripe filling in the true ground state is slightly less 
(0.86 < Ufi < 0.89 hole/Ni depending on the crystal field 
splitting) than the experimental value of one hole per Ni 
ion. This concerns the entire low doping regime, x < 0.3, 
where e ~ x. There may be a few reasons of this discrep- 
ancy. First of all, one has to realize that in a multiband 
model with oxygen orbitals included explicitly the holes 
would be doped primarily to oxygen orbitals, screening 
the local moments at Ni sitesj^i in analogy to Zhang- 
Rice singlets in the cupratesiS An insulating ground 
state could then result from oxygen distortions by the 
Peierls mechanism. Thus, we argue that the present re- 



TABLE III: Comparison of the free energy F (in the units 
of t) for the BC and SC stripe phases in the Eg model, and 
for BC stripe phases in the DDH model, obtained on 64 x 64 
clusters at j: = 1/8 hole doping. The free energy of the most 
stable phase (for BC or SC phases) is given in bold characters 
in each case. Parameters: U — 8t, Jh = 1.5t, E^ = 0, and 
pt = 100. 



model 



DDH 





diagonal stripes 


vertical stripes 


hase 


filled 


half-filled 


filled 


half-fiUed 


BC 


2.5508 


2.5759 


2.5629 


2.5756 


SC 


2.7815 


3.2615 


2.7012 


3.0729 


BC 


2.8172 


2.8135 


2.8395 


2.8354 
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suits could be further improved within a realistic model 
including not only two Cg orbitals with different hopping 
elements, but also orbital polarization at oxygen sites. 
Second, it may be expected that the stripe filling will be 
somewhat changed due to the electron correlation effects 
beyond the HF approximation used here. We believe that 
the present work provides a good starting point for future 
studies of the correlation effects. 

In summary we have analyzed the stripe formation in 
the doped nickelates in the realistic model with degen- 
eracy of eg orbitals, using large clusters. The results 
obtained with this model were compared with the widely 
used doubly degenerate Hubbard model. Even though in 
both models the distance between the DW is inversely 
proportional to the doping, the most stable diagonal 
stripes are found to differ markedly from one model to 
the other one. In the eg model, the stripes are filled and 
nearly insulating, as observed experimentally in a series 
of layered nickelates. In contrast, in the DDH model they 
are half-filled and metallic. These latter stripe phases 



are reminiscent of the stripes observed experimentally in 
largely doped cuprates. These differences have their roots 
in the different structure of intersite hopping terms. As 
the DDH model is closer to t2g than to eg hopping ma- 
trix elements, one might expect that doped insulators 
with t2g orbital degrees of freedom, as in vanadates or in 
ruthenates, would promote different stripe phases than 
those observed in the nickelates. 
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